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ABSTRACT 

Satellite thermal infrared (IR) spectral emissivity data have been shown to be significant for atmospheric research and 
monitoring the Earth’s environment. Long-term and large-scale observations needed for global monitoring and research 
can be supplied by satellite-based remote sensing. Presented here is the global surface IR emissivity data retrieved from 
the last 5 years of Infrared Atmospheric Sounding Interferometer (IASI) measurements observed from the MetOp-A 
satellite. Monthly mean surface properties (i.e., skin temperature T s and emissivity spectra e.) with a spatial resolution of 
0.5x0. 5-degrees latitude-longitude are produced to monitor seasonal and inter-annual variations. We demonstrate that 
surface e, and T s retrieved with IASI measurements can be used to assist in monitoring surface weather and surface 
climate change. Surface e, together with T s from current and future operational satellites can be utilized as a means of 
long-term and large-scale monitoring of Earth’s surface weather environment and associated changes. 
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1. INTRODUCTION 

Earth surface spectral emissivity (e.) is a measure of the relative ability of the surface to emit energy by radiation at the 
frequency v. It is the ratio of energy radiated by its surface to energy radiated by a blackbody at the same temperature. 
The value of Earth surface emissivity depends on the composition of the surface materials and varies from location to 
location and from time to time. Great effort has been devoted toward obtaining and validating Earth’s land surface 
emissivity in a broad wavelength region. Emissivity in the thermal infrared (IR) and microwave regions is particularly 
interesting as it is crucial for achieving other highly-accurate retrieved parameters, the assimilation of radiances in 
numerical weather prediction models, climate modeling/monitoring, and studies of the radiation budget for the Earth 
system. 

Data from satellite remote sensors developed for various environmental sensing applications have been useful for 
deriving surface thermal IR emissivity. For example, emissivity data have been derived using measurements from the 
Moderate Resolution Imaging Spectroradiometer (MODIS), the Advanced Spaceborne Thermal Emission and Reflection 
Radiometer (ASTER), the Spinning Enhanced Visible and Infrared Imager (SEVIRI), the Atmospheric InfraRed Sounder 
(AIRS), and the Infrared Atmospheric Sounding Interferometer (IASI) 1 " 6 . Emissivity derived from satellite 
measurements is an “effective” emissivity of the area where the measurement is taken; it depends on the composition of 
the surface materials in that area. This effective emissivity has been assumed by a linear combination of pure scene 
types 7 s . The vegetation cover and soil moisture of the area can significantly affect the emissivity derived from satellite 
measurements 9 ' 10 . The IR e. derived from IASI measurements has provided spectral information of the surface materials 6 
for investigations carried out and reported herein. The aim of this manuscript is to show global s, climatology from the 
last 5 years of IASI measurements and initial investigation of its seasonal variation over that period of time. The period 
of 5 years is long enough to observe its variation and to indicate and/or recommend that e, and its trend can be used in 
assisting surface weather monitoring and climate studies. 
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2. EMISSIVITY DATA FROM IASI 


The Infrared Atmospheric Sounding Interferometer (IASI) was successfully launched into polar orbit aboard the 
operational polar orbiting MetOp-A satellite on 19 October 2006; significant scientific accomplishments have been made 
with IASI measurements 1112 . IASI is a Michelson interferometer with spectral coverage between 3.6 and 15.5 pm. At 
nadir, the instrument samples data at intervals of 25 km along a cross track, each sample having a maximum diameter of 
~12 km. It possesses a continuous spectral coverage from 645 to 2760 cm" 1 and a high spectral resolution of 0.25 cm" 1 ; in 
other words, its unapodized spectral resolving power (v/8v) is between 2000 to 4000 13 . This higher spectral resolving 
power benefits retrieving e. simultaneously with T s and other atmospheric parameters. An inversion scheme, dealing with 
cloudy as well as cloud-free radiances observed with ultraspectral IR sounders, has been developed and applied to IASI 
data to retrieve simultaneously atmospheric thermodynamic and surface or cloud microphysical parameters 14 . The data 
used in this study are produced with the algorithm described by Zhou et al. 6 . Initial IASI s, evaluations were 
performed 15 ' 16 , providing confidence on the accuracy of the products. 

Here the so-called “monthly mean climate” dataset (hereafter denoted as MMC) is defined as the monthly mean of all the 
years of IASI measurements (i.e., from June 2007 to May 2012), while the “monthly mean” is the monthly mean of a 
certain year (hereafter denoted as MMY). The data have been integrated into a spatial grid of 0.5x0. 5-degrees latitude- 
longitude. IASI s, is spectrally resolved covering all IASI channels. Figure 1 shows the global distributions of December 
nighttime MMC T s and e, at a frequency (v) of 890 cm" 1 (later simplified as £ 89 o). IASI derived e, shows a seasonal 
variation, explained by a combination of surface materials, modified from time to time 6 . In other words, the soil 
moisture, snow/ice cover, and vegetation canopy varies over a different surface backgrounds (e.g., bare soil, sand, and 
rock mixed surface). Moreover, the e, contains its surface material information in different spectral frequencies, identified 
in both laboratory and satellite measurements. For instance, emissivity in the quartz reststrahlen bands (750-1250 cm" 1 ) is 
apparently dominated by the residual mineral quartz 17,18 . Figure 2 illustrates a few £, representing Earth surface materials, 
which are groups of some similar materials measured in the laboratory 17 . For example, ice £, is smaller than that of water 
in IR region. Although sea ice s, can vary depending on sea ice type, the £, difference between sea ice and seawater is 
small but significant, especially in the shortwave region. For instance, using the MMC data of December, the £2400 is 
0.981 and T s is 296.76° K at a location of 0.0-0. 5° N Latitude and 149.5-150.0° W Longitude (over seawater); while the 
£2400 is 0.955 and T s is 266.74° K at a different location of 59.5-60.0° S Latitude and 149.5-150.0° W Longitude (most 
likely over sea ice). This helps us to understand why the s, shown in Figure lb is slightly smaller in the higher latitude 
regions over the water where sea ice is present. 
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Figure 1. IASI monthly mean climatology: December nighttime (a) surface skin temperature (T s ) and (b) 

emissivity at 890 cm" 1 (£ 89 o). 


3. EMISSIVITY SEASONAL VARIATIONS AND MONITORING 

Because the nature of emissivity does not change rapidly, monthly mean emissivity is used and presented for analysis 
and monitoring. Emissivity seasonal variation has been indicated in the satellite measurements. For illustration, Figure 


3a plots the MMC e, measured from Utah over the Rocky Mountains (40.0-40.5° N Latitude; 109.5-110.0° W Longitude; 
~1.8 km elevation). Channel emissivity Shoo and e 2 4 oo at 1100 cm" 1 and 2400 cm' 1 , respectively, for each month is plotted 
in Figure 3b, showing that the higher emissivity appears in February and July while the lowest one is in April. Surface T s 
is plotted in Figure 3c. It implies that snow/ice covers the surface and mostly melts down by April when the emissivity 
reaches the lowest point with its associated T s of 276.64° K and fresh-vegetation cover is still at a minimum in the early 
spring. It increases from April to July when soil moisture might be high and fresh-vegetation cover increases. It then 
decreases from August to November as vegetation cover decreases. 



Figure 2. Samples of emissivity spectra representing Earth surface materials. 
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Figure 3. (a) Monthly mean s, climatology measured from 40.0-40.5° N Latitude and 109.5-1 10.0° W 
Longitude, (b) Channel emissivity at 1100 and 2400 cm" 1 as a function of time, (c) Surface skin 

temperature as a function of time. 

Although both winter (February) and summer (July) emissivity is relatively high and their difference is small, it is 
noticed that their spectral shapes are different. For instance, their difference depends on the frequency (or wavenumber). 
For the case of spring to summer, we assume that flash tree leaves (e.g., Quaking Aspen, common in the area) and some 
other native vegetation start to grow in the spring. Vegetation cover reaches a maximum by July. We use laboratory- 
measured white pine tree emissivity, roughly representing a part of the surface materials, and April-measured emissivity 
to reproduce a combined emissivity to fit July-measured emissivity. Figure 4a shows the combination of 68% of April 
IASI-measured emissivity and 32% white pine emissivity can closely match July IASI-measured emissivity. Similar 
analysis has been done for the winter-to-spring case assuming that mixed ice-snow partially covers the ground during the 
winter season. The April IASI-measured and laboratory-measured ice emissivity are used to make a combined emissivity 


to match February IASI-measured emissivity. Figure 4b shows the combination of 54% of April IASI-emissivity and 
46% ice emissivity can match February IASI-emissivity. This implies that the emissivity in this region can be used to 
monitor a relative amount of ice/snow during the winter season and a relative amount of native vegetation during the 
summer season. 




Figure 4. (a) Vegetation increasing from April to July; July eis roughly reproduced by a linear 
combination of April e, and white pine e„ (b) Snow/ice covers decreasing from February to April; February 
s. is roughly reproduced by a linear combination of April e. and ice 8.. 

We chose IASI channel emissivity at the wavenumber of 975 cm" 1 (8975) for a trend analysis and used a few locations to 
demonstrate that its variation behavior is location dependent. Figure 5 plots monthly mean (a) T s , (b) T s anomaly ( 4 > T s), 
(c) 8975, and (d) 8975 anomaly (4> e 97 5 ) from June 2007 to May 2012. The 4 > Ts (or (^975) is defined as the monthly mean of the 
year (MMY) deviated from the monthly mean climatology (MMC). The first location is from the Utah Rocky Mountains 
as described above. Figure 5c shows the monthly mean 8975 (filled red triangle) as a function of time plotted with 
monthly mean climatology (filled blue circle). The difference (MMY-MMC) or anomaly 4)5,975 is plotted in Figure 5d. It 
is clearly evident that “snow/ice” cover varies; for example, January to March 2012 had higher T s and lower S975 (less 
snow/ice) than those in the same months during previous years. More snowfall during winter of 2010-2011 led to extra 
soil moisture and native vegetation during the following 2011 spring and summer. In the figure, the standard deviation 
errors (shown with error bars) are estimated from a single field-of-view (SFOV) retrieval 6 and then divided by the square 
root of the number of SFOV retrievals used for a monthly mean in the spatial grid of 0.5x0. 5-degrees latitude-longitude. 
The absolute deviation (i.e., bias) is assumed to be similar in both MMC and MMY data; therefore, the biases roughly 
cancel each other in their differences (i.e., anomaly (|)). The red line plots a linear fit of the anomalies, indicating the 
anomaly trend during this period. At this point, this trend may not be useful in drawing any tendency conclusions, but it 
certainly can be used for monitoring its variations as more data will be added on from upcoming measurements. 

It is worth pointing out that s, variation may not be observed (i.e., no variation) in some places. For example, the monthly 
mean 8, over water is not changed until the water is cold enough to form ice. In Figure 6, the left column is from the 
equatorial Pacific (0.0-0. 5° N Latitude; 139.5-140.0° W Longitude); Figures 6al and 6b 1 show that T s is almost 
seasonally independent except during the period when La Nina and El Nino are currently evident. The (Jit, changed from 
warm to cool between April to May 2010 (Figure 6b 1 ), implying that the tropical Pacific switched from El Nino to La 
Nina. This was also tracked and confirmed by the NASA/European Ocean Surface Topography Mission/Jason-2 
oceanography satellite. As El Nino plays a significant role in tropical drought occurrence 19 , a sample from Australia is 
chosen (right-column of Figure 6: 24.5-25.0° S Latitude; 145.0-145.5° E Longitude) to illustrate e, variation associated 
with droughts during El Nino years and greater rainfalls during La Nina events. Figures 6c2 and 6d2 plot 8975 and ^975, 
respectively. The 8 increases when soil moisture increases as greater rainfalls occurred; in contrast, the s decreases when 
soil moisture decreases as droughts occurred 10 ' 19 . 
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Figure 5. Sample from 40.0-40.5° N Latitude and 109.5-1 10.0° W Longitude: Illustration of monthly 
mean (filled red triangle) (a) T s and (c) 8975 against their monthly mean climatology (filled blue circle); and 
(b) T s anomalies and (d) 8975 anomalies plotted to observe their trends (see text). 
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Figure 6 . Same as Figure 5 but at 2 other locations: The El Nino (or La Nina) phenomena shown in (bl) 
causes the land emissivity change shown in (d 2 ) related to the soil moisture, associated with droughts (or 

great rainfalls). 


4. CONCLUSION 

We have developed a global climatology atlas of IR e, and T s based on the last 5-years of IASI measurements. The IR e. 
retrieved with IASI measurement is sensitive to surface material. Different surface material is not only reflected in an 
absolute channcl-emissivity value at a certain frequency but also in its spectral shape. Retrieved s. shows a temporal 
variation as its surface and/or surface weather changes. Relatively speaking, the change of surface cover (e.g., 
vegetation, snow/ice) and soil moisture can be detected with its e. depending on the surface background. Month-to-month 
a variation is clearly evident and it can be significant in some areas such as deserts and bare land. Temporal variation 
indicates a seasonal diversity; El Nino and/or La Nina effects are evident not only on the water (i.e., equatorial Pacific 


T s ), but also on the land surface (i.e., Australia land e,). The different e. for sea ice and seawater might be used, in 
conjunction with its associated T s , to monitor the change of global sea ice amounts. Initial findings of emissivity seasonal 
and inter-annual variation phenomenon at a few locations are reported herein. As the retrieval algorithm continues to be 
improved to provide higher-quality retrieval data, more investigation and studies with the data during a longer time 
period are strongly desired to draw explicit conclusions. However, a long-term monitoring of the Earth’s surface 
environment can be critical in assisting weather and climate research. Satellite-measured surface emissivity can be used 
as a means of long-term and large-scale monitoring of the Earth’s environment for climate studies. 
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